Wavelength calibration for OSIRIS/GTC* tunable filters 

J. Mendez-Abreu^'^, J. Sanchez Almeida^'^, C Muiioz-Tuiion^'^, J. M. 
Rodriguez-Espinosa^'^, J. A. L. Aguerri^'^, D. Rosa Gonzalez ^, Y. D. Mayya ^, O. Vega ^, 
O ■ R. Terlevich ^'^, E. Terlevich^'^, E. Bertone ^, L. H. Rodriguez- Merino ^ 



O 
(N 



ABSTRACT 



OSIRIS (Optical System for Imaging and low Resolution Integrated Spec- 
troscopy) is the first light instrument of the Gran Telescopio Canarias (GTC). 
It provides a flexible and competitive tunable filter (TF). Since it is based on 
(-H ' a Fabry-Perot interferometer working in collimated beam, the TF transmission 

Ph' wavelength depends on the position of the target with respect to the optical axis. 

O ■ This effect is non-negligible and must be accounted for in the data reduction. 

H ■ 

'^ ' Our paper establishes a wavelength calibration for OSIRIS TF with the accu- 

^ ■ racy required for spectrophotometric measurements using the full field of view 

(FOV) of the instrument. The variation of the transmission wavelength A(-R) 

> : across the FOV is well described by X{R) = A(0)/'yiTcRZ^. where A(0) is 

O '• the central wavelength, R represents the physical distance from the optical axis, 

^ ! and /2 = 185.70 ± 0.17mm is the effective focal length of the camera lens. This 

1^ [ new empirical calibration yields an accuracy better than 1 A across the entire 

O I OSIRIS FOV (~8'x8'), provided that the position of the optical axis is known 

within 45 /xm (= 1.5 binned pixels). We suggest a calibration protocol to grant 
such precision over long periods, upon re-alignment of OSIRIS optics, and in dif- 
k/< I ferent wavelength ranges. This calibration differs from the calibration in OSIRIS 

H I manual which, nonetheless, provides an accuracy < lA for R < 2'. 

Subject headings: Data Analysis and Techniques — Astronomical Techniques 



This work is based on observations made with the GTC operated on the island of La Palma by Grantecan 
in the Spanish Observatorio del Roque de los Muchachos. 
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Introduction 



Tunable filter (TF) instruments are becoming standard tools for mapping physical prop- 
erties in extended astronomical sources. Most modern large telescopes have TFs based on 
Fabry- Perot interferometers (FPI; or etalons) such as the 11m Southern African Large Tele- 
scope (SALT; Rangwala et al. 2008) or the Magellan-Baade 6.5 m telescope (Veilleux et 
al. 2010). The tunable imaging era began with the instrument of Atherton & Reay (1981) 
and many others have further developed the technique (Bland & Tully 1989; Brown et al. 
1994). The current instrumental design is mainly based on the Taurus Tunable Filter (Bland- 
Hawthorn & Jones 1998) formerly mounted at the Anglo- Australian Telescope and afterward 
moved to the William Herschel Telescope at La Palma. The power of these instruments to 
perform spectrophotometric studies is clear. They allow for spectrophotometry with an un- 
limited number of wavelengths that can be scanned within a wavelength range. These etalon 
based TFs provide a central monochromatic field (MF). Outside the MF, and due to the use 
of etalons in collimated beam, there is a significant shift of wavelength that depends on the 
distance to the center. The instruments are designed to provide the largest possible MF, the 
wider possible wavelength range, and the narrower possible transmission band-pass. Besides, 
a proper calibration corrects the wavelength displacement with radius, thus recovering much 
of the field outside the MF to be used as an extended MF. 

Among the large facilities hosting such instruments is the 10 m Gran Telescopio Canarias 
(GTC), located at the Roque de los Muchachos Observatory, that started its operations in 
the first semester of 2009 with OSIRIS (Cepa et al. 2005; Cepa 2010) as the first light 
instrument. OSIRIS stands for "Optical System for Imaging and low Resolution Integrated 
Spectroscopy". Besides conventional imaging and spectroscopy, it provides the additional 
capability of a TF observing mode working in the visible. It covers the wavelength range 
0.365-1.0 /xm with a nominal FOV of 7.8'x8.5'. Further commissioning tests establish the real 
un- vignetted FOV as 7.8' x 7.8' (see OSIRIS commissioning webpage^). OSIRIS specifications 
are extensively described in the manual provided by the instrument team (OSIRIS TF user 
manual^ - this document is referred to along the paper as OSIRIS manual). In particular, 
it provides the wavelength calibration procedure for radii outside the MF. 

The OSIRIS FPI is an ideal instrument to map emission line regions in galaxies. It 
provides a high throughput with narrow transmission bandpass. The central wavelength 
is adjustable to account for the galaxy recession velocities, and the field of view (FOV) 



^ http://www.gtc.iac.es/en/pages/instrunicntation/osiris/data-coiiimissioning.php 

^http://www.gtc.iac.es/en/pages/instrunicntation/osiris/osiris2.php^Useful_Docuincnts; current version 
dated June 28, 2009 
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is reasonably large with a small pixel scale (0'.'125). These properties together with the 
collecting power of a 10 m telescope, designed to get the best possible natural image quality, 
prompted us to design the Local Universe Survey (LUS^), a research program optimized for 
OSIRIS. LUS aims at studying the star-formation history of a complete sample of nearby 
galaxies. We planned to obtain high signal-to-noise, high angular resolution, narrow band 
maps of all galaxies inside a volume of 3.5 Mpc radius. All irregulars and spirals inside 
a volume of 11 Mpc were also included, together with a sample of Virgo cluster galaxies. 
These observations suffice to produce a detailed unique description of the gaseous and stellar 
content of nearby galaxies. LUS has been defined within the ESTALLIDOS^ project in close 
collaboration with the Instituto Nacional de Astrofisica, Optica y Electronica in Mexico. 
Significant preparatory work was carried out in order to prove the feasibility of LUS; using 
the OSIRIS features described in the manual, we made simulations to estimate the exposure 
times, to choose precisely the filters (broad band and TF), and to optimize the observing 
procedure. Before submitting LUS as an ESO-GTC large program, we decided to test our 
simulations and procedures with real data. 

The object selected to calibrate the photometric accuracy was MlOl, an extended and 
well-known nearby galaxy (7.38 Mpc, Rizzi et al. 2007), that exceeds the size of OSIRIS 
FOV. Observations were carried out to cover a large region of the southwestern part of 
the galaxy, containing also the well-known giant extragalactic H II region NGC 5447. The 
central wavelengths of the TF were chosen to map emission lines such as Ha or the sulphur 
doublet ([S II]AA6717, 6731), which allow to derive star-formation rates (SFRs) and electron 
densities of the star-forming regions. Even if not optimal, observations were made using the 
narrowest bandpass of FWHM 18 A available at the time. The precise central wavelength 
of the filter needed to separate [N II] A6584 from Ha was also tested, as distinguishing the 
two lines is critical to be able to use diagnostic diagrams like BPT (Baldwin et al. 1981). 
The basic data reduction was carried out using our own pipeline. Then, we reconstructed 
monochromatic emission line images taking into account the wavelength calibration provided 
in OSIRIS documentation. From this analysis, we found that Ha fiuxes and SFRs were 
globally consistent with values from the literature. However, the emission line ratios for 
both [NII]A6584 over Ha and the sulphur doublet ([S II]AA6717, 6731) were far from any 
observed value and whatever H II region model. In an attempt to identify the problem, we 
started from scratch by using spectral calibration lamps to test wavelength shifts beyond 
the MF. The result of our calibration is the content of this paper. We have figured out that 
the prescriptions given in the manual suffice for the MF, while the new calibration in this 



'^http://www.inaocp.nix/ gtc— lus/ 

'* http://estallidos.iac.es/estallidos/Estallidos.jsp 
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paper allows for using OSIRIS in its full capability. It permits spectrophotometry in the 
full FOV with a wavelength accuracy as good as that in the MF. In our view OSIRIS is an 
extremely powerful instrument which could be seriously limited if only the MF can be used 
for spectrophotometry. Although our work was designed to recalibrate our fields and being 
able to recover reliable emission line fluxes over the entire OSIRIS FOV, it is of interest 
beyond the original scope, and may broaden the community of OSIRIS users. 

The paper is organized as follows: in § 2 we introduce the basic optical concepts required 
to understand both the nature of the problem and the calibration procedure. § 3 shows how 
the original wavelength calibration needs refinement to provide an accuracy better than lA 
in the full FOV. Based on the data introduced in § 4, the new calibration is worked out in § 5. 
In § 6 we recommend a wavelength calibration protocol that goes all the way from a detailed 
test of the long-term stability of the system, to a black-box recipe for urgent calibration. 
The conclusions of the work are summarized in 5 7. 



2. Basics of OSIRIS/GTC tunable filter optics 
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Fig. 1. — Schematic of the OSIRIS/GTC optical layout. A light beam, in red, crosses 
the system from left to right, finding the GTC telescope, the collimator, the Fabry- Perot 
interferometer, the camera, and OSIRIS focal plane. The different symbols are defined in 
the main text. 



In essence, the optical layout of OSIRIS/GTC TFs is made of three lenses and a FPI, 
as sketched in Fig. 1. The first element, of effective focal length /o, represents the telescope 
(GTC). Then a lens, of effective focal length /i, collimates the light coming from the GTC 
focal plane, so that the FPI is illuminated in collimated beam. As represented in Fig. 1, 
a point in the focal plane of GTC, at a distance D from the optical axis, enters the FPI 
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interferometer with an angle 9. Then 

D = fitanO = /otana ~ foa, (1) 

where a represents the angular distance to the optical axis of the source on the sky. The 
wavelength transmitted by a FPI depends on the incidence angle as 

A(^) = A(0)cos^, (2) 

where A(0) is the central wavelength (e.g., Born & Wolf 1980; Beckers 1998; Veilleux et al. 
2010). The final optical element, the camera, has an effective focal length /2, so that, 

R = f2tan9, (3) 

where R is the position of the source on OSIRIS focal plane. Combining Eqs. (2) and (3), 
one gets the TF wavelength calibration equation in terms of the physical distance to the 
optical axis R. It is given by 

Mu) ^ ^i£L=, (4) 

or equivalently, employing Eq. (1), it can be expressed in terms of the angular distance 
projected on the sky 

w . A(0) 

A(") = / =• (5) 

l+(af 
If a is small enough, then Eq. (5) can be expanded to first order giving rise to 

A(«) ~ A(0)(l-Ca2), (6) 

n — ■^o 

and equivalently in terms of the physical distance 

\{R) - A(0)(l-C*i?2), (7) 

1 



C* 



2/1 • 



The plate scale on OSIRIS focal plane relates the effective focal lengths of the three 
optical elements. Combining Eqs. (1) and (3), the plate scale S turns out to be 

da_ /i ^^^ 



dR kf2 



-6- 

which we use to link the two calibration constants, 

C* = CS^. (9) 

OSIRIS manual employs Eq. (6) for wavelength calibration, therefore, it assumes the 
FOV to be small enough for the first order expansion to hold. In addition, the calibration is 
performed in angular coordinates on the sky, so that uncertainties in the plate scale affect the 
wavelength calibration. We propose using Eq. (4) instead. It makes the calibration sensitive 
only to the uncertainties of the focal length of the camera /2, and it considers all the terms 
of the expansion. As we show in the next section, high order terms are needed to calibrate 
the field outside the MF. 

3. Motivation for an updated calibration 

Inconsistencies in emission line ratios prompted us to reconsider OSIRIS wavelength cal- 
ibration (§ 1). The importance of this effect can be disclosed and quantified using monochro- 
matic images of known wavelength, as those provided by the GTC wavelength calibration 
unit. When the instrument is illuminated with a monochromatic beam of known wavelength 
Aline, the ring appearing in the focal plane traces where OSIRIS transmission matches Aune 
(see Fig. 2). In principle, the position of the optical axis and the tuned central wavelength 
A(0) are given, therefore, one can easily measure the difference between the true wavelength 
Aline, and that provided by the nominal calibration in OSIRIS manual (Eq. (6)). If the ring 
is observed to have a radius R (Fig. 2), then the error in wavelength is given by 

AA = Aline - KR) = Aline - A(0) (1 -CS^R"), (10) 

where C = 7.9520 x 10~^arcmin~^ is the calibration constant appearing in the manual, and 
S is the nominal plate scale (0'.'125 pixel"^; see Table 1 with OSIRIS properties). 

The black dots in Fig. 3 show A A for a suite of 17 monochromatic images used for 
calibration (details are given in § 4). Note that the error in the true wavelength can be as 
large as 6 A at the FOV outskirts {R ~ 4'; see Fig. 3). If rather than using the approximate 
Eq. (6), theoretical wavelengths are computed using the full Eq. (5), then one obtains the 
triangles in Fig. 3. The errors have been somewhat reduced, but they still remain as large 
as 4 A at 4'. Two main conclusions can be drawn from Fig. 3. First, the original calibration 
procedure does not provide correct wavelengths within 1 A when R > 2'. Second, in order 
to achieve such an accuracy, one has to use the full Eq. (4) combined with a new calibration 
constant more precise than that inferred from the nominal values of the focals. 



Table 1. Summary of the basic optical parameters for OSIRIS/GTC wavelength 

calibration. 



Physical Parameter 



Symbol 



Nominal Value 



Best Estimate 



GTC eff. focal length 
Collimator eff. focal length 
Camera eff. focal length 
Calibration constant 
Calibration constant 

OSIRIS pixel size (1x1) 
OSIRIS pixel size (2x2) 
OSIRIS plate scale (1x1) 
OSIRIS plate scale (2x2) 
OSIRIS/GTC eff. focal length 
OSIRIS center position (2x2) 



/o 


170 m [a] 




165.3 ± 0.2 m [h] 


/l 


1.240 m [b] 




1.275 ± 0.001m [i] 


/2 


181mm [b, c] 




185.70 ± 0.17 mm 


C = fi/{2fl) 


7.9520 X 10"**arcmin-2 


[d] 


(7.517 ± 0.014) X 10-''arcmin-2[j] 


C*= 1/(2/1) 


15.262 m-2 




14.499 ± 0.027 m-2 




— 




(3.2624 ± 0.0057) x 10"^ pixcr^ 




15 /im [a, b, c, d] 




— 




30 /im 




— 


S=/l/(/0/2) 


0'.' 125 pixel- 1 [a, b,d,e] 




0'.'127pixer^ [k] 




0'.'250 pixel- 1 




— 


/0/2//1 


24.752 m [b, g] 




— 


(xci, yd) 


(1059± 2, 983± 2) [f] 




(1053.6 ± 2.2,979.5 ± 1.8) 


(Xc2, yc2) 


(-1± 2, 981± 2) [f] 




(-9.8 ± 2.7,975.8 ± 1.1) 



''Cobos et al. (2002). 

''OSIRIS project wcbpage http://www.iac.es/project/OSIRIS/. 

'^Cepa et al. (2005). 

''Value included in OSIRIS/GTC TF user manual, and derived from the quoted nominal values for /o and /i . 

"Cepa (2010). 

^OSIRIS/GTC TF webpage http : //www. gtc . iac . es/en/pages/instrumentation/oslris .php#Tunable_Filters . 

sit does not agree with the other nominal values - it is obtained for /o = 169.57 m. 

'^ Value calculated using Eq. (8) assuming the nominal values for S and /i, and the best estimate for /2. 

'Value calculated using Eq. (8) assuming the nominal values for S and /o, and the best estimate for /2. 

J Value calculated assuming the S =0''125pixcl--'^ nominal value. For other plate scales, see Eq. (14). 

"^ Obtained from astrometry and in agreement with laboratory tests. 
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4. Data 

Figure 3 and the calibration worked out in § 5 are based on a sample of 17 monochromatic 
images produced by arc lamps from the Instrument Calibration Module of GTC (ICM). They 
were obtained in two different runs and using the OSIRIS red arm. Their main properties 
are summarized in Table 2. The images were taken by the GTC staff who kindly provided 
them upon request. The OSIRIS FOV is imaged onto two CCDs that are separated by a gap. 
The distance between the CCDs is approximately 9'.'2. The two OSIRIS science CCDs are of 
type E2V 44-82, back illuminated, with 2048 x 4096 pixels each when using a 1x1 binning. 
OSIRIS was tuned to different central wavelengths in order to produce the characteristic 
ring for every monochromatic line (Fig. 2). They were obtained using the standard mode 
of operation with a binning 2x2. The separation between any two consecutive orders of 
interference is called the free spectral range. In order to select the desired transmission 
order of the FP a blocking filter is required. The filters used in this study are also shown in 
Table 2. 

In addition to the above monochromatic exposures, three images with strong sky emis- 
sion rings were also requested in order to compare the wavelength calibration obtained from 
the ICM and from the sky. They are also listed in Table 2. 



5. Wavelength calibration 

As we argue in § 3, improving the original calibration demands both including all the 
terms neglected in the approximate Eq. (6) and determining the calibration constants directly 
from observations. We measured C* because it does not require knowing the plate scale S", 
therefore, we make the calibration insensitive to errors in S. Besides, C follows directly from 
C* and S via Eq. (9). 

Calibrating an image is assigning a wavelength to each point of the FOV using Eq. (4). 
Three different unknowns must be set, namely, the position of the center, origin of distances, 
the calibration constant, and the central wavelength. This section describes a precise determi- 
nation of the two first ingredients using the data described in § 4. The central wavelength is 
determined within 1 A as part of the regular OSIRIS tuning protocol. This claim is confirmed 
in § 5.2 since it also affects the empirical determination of C* . The three unknowns have un- 
certainties that determine the calibration accuracy. The propagation of these uncertainties 
into the final wavelength is studied in § 5.3. 
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Table 2: Characteristics of the cahbration data. 



A(0) 


-^line 


Species 


TF FWHM 


Blocking Filter 


Observing 


(A) 


(A) 




(A) 


(nm/nm) 


Date 


6598.9 


6532.9 


Nel 


18 


f657/35 


10-02-2011 


6630.0 


6598.9 


Nel 


18 


f657/35 


10-02-2011 


6650.0 


6598.9 


Nel 


18 


f657/35 


10-02-2011 


6950.0 


6929.5 


Nel 


16 


f694/44 


03-03-2011 


6975.0 


6929.5 


Nel 


16 


f694/44 


03-03-2011 


7000.0 


6929.5 


Nel 


16 


f694/44 


03-03-2011 


7020.0 


6929.5 


Nel 


16 


f694/44 


03-03-2011 


7050.0 


7032.4 


Nel 


16 


f694/44 


03-03-2011 


7070.0 


7032.4 


Nel 


16 


f709/45 


03-03-2011 


7090.0 


7032.4 


Nel 


16 


f709/45 


03-03-2011 


7100.0 


7032.4 


Nel 


16 


f709/45 


03-03-2011 


7120.0 


7032.4 


Nel 


16 


f709/45 


03-03-2011 


7650.0 


7635.1 


Arl 


18 


f770/50 


03-03-2011 


7680.0 


7635.1 


Arl 


18 


f770/50 


03-03-2011 


7700.0 


7635.1 


Arl 


18 


f770/50 


03-03-2011 


7720.0 


7635.1 


Arl 


18 


f770/50 


03-03-2011 


7740.0 


7635.1 


Arl 


18 


f770/50 


03-03-2011 


6900.0 


6834.4-6841.9 


telluric OH 


12 


£666/36 


12-02-2011 


6925.0 


6871.0-6880.9 


telluric OH 


12 


f666/36 


12-02-2011 


6950.0 


6864.0-6871.0 


telluric OH 


12 


f680/36 


12-02-2011 
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Fig. 2. — Image of a Nel spectral calibration lamp illuminating the OSIRIS detectors. It 
forms a ring which we have outlined with a green circle of radius R centered at X. OSIRIS 
has been tuned to a central wavelength of 7070 A, whereas the Nel emission occurs at 7032.4 
A. Binned coordinates (1,1) and (1049,2051) are also drawn on the image. 

5.1. TF optical center 

The nominal TF optical center is located in the pixel (2118, 1966) from the pixel (1,1) 
of CCDl, including the 50 pixels of overscan. Thus, assuming the nominal gap between the 
CCDs of 72 pixels^, the center lies within the gap of the two CCDs (see Fig. 2). As pointed 
out before, our images have been taken with the standard 2x2 binning. Therefore, we might 
expect a center positioned at binned pixels (1059,983) - see Table 1. Taking into account 
that the distance between the CCDs is ~36 binned pixels, the TF optical center should be 
placed 1 pixel to the left of column 1 of CCD2. 



^http://www.gtc.iac.es/en/pages/instrunientation/osiris.php:^Tunable_Filters. 
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Fig. 3. — Differences between the true wavelengths and that derived using the nominal 
OSIRIS calibration procedure. Each point has been inferred from one image similar to that 
in Fig. 2. The discrepancies increase with distance from the optical axis R (= Radius). The 
black dots correspond to the nominal calibration based on the approximate Eq. (6), whereas 
the blue triangles used the full Eq. (5). Abscissas are given in pixels (bottom) and arcsec 
(top), the latter calculated using the nominal plate scale. Vertical dotted lines show the 
radius of the MF (~2') and the non- vignetted FOV (~4'). 



We have measured the optical center for our calibration lamps using the following proce- 
dure: first we generate a grid of possible TF optical pixel centers (xtriai,Z/triai) independently 
for both CCDs. Since the gap between CCDs is not constant (it varies from top to bottom 
due to both a rotation and a 2 pixels vertical shift of CCD2 with respect to CCDl), the solu- 
tion of assuming independent centers bypasses the difficulty of knowing the relative position 
of the two CCDs. Then, for every trial position on each CCD, the calibration lamp image 
was divided in 8 regions of 45° (4 in each CCD). Each region is azimuthally averaged and the 
distance of the ring to the assumed center is measured. Due to both possible irregularities 
in the optical system and the offset between the true and trial centers, the measurements 
of the ring radii for each region are not the same. We assign the best center from our grid 
of trial values looking for the value that, providing the same mean center within the errors 
for both CCDs, minimizes the scatter among the radii of the eight different sectors. If the 
monochromatic images were noiseless and the ring perfectly circular, our center-finding algo- 
rithm would provide the correct value. In practice, however, the procedure has some errors 
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that we estimate by comparing centers obtained from different rings. 

Figure 4 shows the best centers obtained for each cahbration lamp image, and the mean 
value over our whole sample. They turn out to be 



(a^ci, 1/ci) = (1053.6 ± 2.2, 979.5 ± 1.8), 

for CCDl, and 

(a:c2, yc2) = (-9.8 ± 2.7, 975.8 ±1.1), (11) 

for CCD2, where each center refers to its own CCD. The error bars in Eq. (11) represent the 
standard deviation among the values obtained from the different rings, therefore, assuming 
the errors of the different measurements to be independent, the error bars of the average 
values should be scaled down by a factor of the order of the square root of the number of 
rings (~4, e.g., Martin 1971). 

In the following we use the mean values as the centers for every calibration lamp. The 
nominal values of the centers are also included in Fig. 4 as triangles. The nominal center 
of CCDl may be marginally in agreement with our empirical determination. However, the 
center of CCD2 is clearly inconsistent with its nominal value. 



5.2. Calibration constants 

A trivial manipulation of Eq. (4) renders 

[||]^^lH-/.-fi^ (12) 

thus, a straightforward linear regression allows us to derive the value of the effective focal 
length of the camera lens (/2) given the positions on the image (i?) of lines of known wave- 
length (A(i?) = Aline)- Figure 5 shows the data and fit leading to the empirical determination 
of /2. We have measured the distance from the center of the system (obtained in § 5.1) to the 
positions of the rings produced by the 17 different monochromatic images described in § 4. 
As it was done in the previous section, this distance was measured 8 times corresponding to 
the 8 different sectors in which each ring was divided. They are plotted in Fig. 5, showing 
a very small scatter that provides confidence on the procedure. The best value for /2 was 
derived by fitting a straight line to the points in Fig. 5 using a standard linear least squares 
minimization routine. Once the slope is properly transformed assuming 1 binned pixel = 
30 /xm (Table 1), the best fit gives 
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Fig. 4. — OSIRIS optical centers obtained using our sample of calibration lamps (black dots). 
The two CCDs are shown separately (CCDl, top; CCD2, bottom). The blue stars represent 
the average value of the centers, which we take as reference for all lamps. The red triangles 
show the nominal center positions. The limits of the plots correspond to the parameter space 
covered by the trial values used to determine the position of the centers. Positions are given 
in binned pixels on the chips. Abscissas larger than 1049 or negative imply centers outside 
the CCDs. 



/a = 185.70 ± 0.17 mm. 



(13) 



The fit also provides the y— intercept. It turns out to be consistent with one (0.99995 ± 
0.00006) supporting the use of Eq. (12) to represent the wavelength variation along OSIRIS 
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FOV. 

Note how the empirical /2 differs from the nominal value of 181 mm (Table 1). This 
difference in the effective focal length of the camera lens was not completely unexpected. 
From astrometry and laboratory tests we know that the actual plate scale of OSIRIS is 
~0'.'127 pixel"^ rather than the nominal (y.'125 pixel"^ (Table 1). This difference can only be 
accounted for if the nominal focal lengths differ from the actual values (see Eq. (8)). 
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4x'0^ 6x10^ 



1 xlO^ 
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Fig. 5. — Fit used to derive empirically the value of the camera lens effective focal length (/2). 
Black crosses represent the values measured using the calibration lamps, one for each region 
of every ring image. The blue line represents the best fit to our measurements. Hatched 
areas indicate the range of possible values for the three sky rings. The black dotted line 
shows the result to be expected using the nominal value for /2, which is way off the observed 
values. 



The hatched areas in Fig. 5 show the position of three rings from sky lines also measured 
in this work (§ 4). The measurements are not represented as points but as regions, since 
they are not single lines but series covering a fairly broad spectral range (Table 2). We have 
considered the wavelengths of the sky lines covered by each filter as an error in the ordinate 
axis. In the abscissa axis, the errors correspond to the position of the 8 different sectors 
in the image. It is worth pointing out that even if the chosen sky lines are not suitable to 
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be used for wavelength calibration, their position on OSIRIS focal plane are fully consistent 
with the new calibration carried out in this work (Fig. 5). 

Once /2 has been obtained, and in order to compare with the nominal values, we have 
derived the calibration constants C* and C using Eqs. (7) and (9), 



C* = 14.499 ± 0.027 m^^ = (3.2624 ± 0.0057) x lO^^pixel'^ 
C = (7.517 ± 0.014) X 10~^arcmin-2 [S/(y!l25pixer^]-^. (14) 

C uses S = 0'.'125 pixel"^ for the sake of comparison with the nominal values, although the 
wavelengths calibrated using C in Eq. (5) are independent of the plate scale. The errors in 
C and C* have been calculated by propagating the error of /2 (Eq. (13)). In the same way, 
the ratio between the GTC and OSIRIS collimator effective focal lengths can be determined 
yielding /i//o = (7.5024 ± 0.0070) x lO'^ with S = 0'.'125 pixel'^ 
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Fig. 6. — Wavelength difference in A between the actual wavelength of the calibration lines 
and that predicted by the new calibration. Black dots represent the measurement for each 
ring sector of every calibration lamp. The black solid lines show the uncertainty caused by 
the error in /2. The red dotted line indicates the deviation from the actual wavelength of 
the calibration lines given by nominal calibration. The scale on top has been calculated 
using the nominal plate scale. Vertical dotted lines show the radii of the MF (~2') and the 
non- vignetted FOV (~4') 
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Figure 6 represents the difference between the true wavelengths of the hnes and that 
inferred using our cahbration. It is equivalent to Fig. 3 but using the new calibration. 
The reason for including this figure is threefold. First, it shows how the systematic trend 
present in Fig. 3 has disappeared, meaning that the new calibration holds throughout the 
whole OSIRIS FOV. Second, the observed scatter provides an upper limit to the error in 
tuning the central wavelength of OSIRIS. As one can see from Table 2, all but one of the 
monochromatic images were taken with a different A(0). Any random error existing in setting 
A(0) would be transferred to A A as a random fluctuation of amplitude similar to that of A(0) 
(see Eq. (10)). The scatter in Fig. 6 does not exceed lA, and so does the setting of A(0). 
Finally, Fig. 6 also shows A A considering the original calibration (the red dotted line). The 
error remains within the acceptable limit of AA < 1 A when R < 2', which approximately 
coincides with OSIRIS MF. 



5.3. Error budget 

Three main sources of uncertainty burden the wavelength calibration, namely, errors in 
the central wavelength (5A(0), errors in the calibration constant 6C*, and errors in the center 
position 6xc, and 6yc. The error in the central wavelength has not been explicitly treated 
before but, according to the arguments given in the previous section, it is within the nominal 
lA set by the procedure to tune OSIRIS (see OSIRIS manual). 

Using Eq. (4), it is possible to propagate these errors onto the wavelength calibration 
6X{R). The result of this exercise turns out to be 

X{R) A(0) ' ^ ^ 

for the central wavelength error (5A(0), 



5\{R) -5C* 



X{R) 2C* 

for the calibration constant error 6C*, and finally, 

6X{R) (5a;c cos 
X{R) ^ R~ 



A^(0) 



A2(0) 



(16) 



;i7) 



for the error of the center in the x-axis {6xc). The variable stands for the azimuth of 
the point relative to the center, and it has been introduced for convenience since using 
I cos 01 = 1 in Eq. (17) grants an upper limit to SX. The error due to uncertainties in the 
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vertical direction is formally identical to that for the horizontal direction replacing 5xc cos 
with 5yc sin 0. 

In order to illustrate the relative importance of the various sources of error, Fig. 7 shows 
the typical 5\{R) to be expected with our calibration, assuming the central wavelength to 
be tuned to Ha (6563 A). The dotted line corresponds to the effect of lA uncertainty in the 
determination of the central wavelength. It sets a minimum threshold to the calibration error 
in the full OSIRIS FOV. The dashed line considers 5C* = 0.027 m~^, which corresponds to 
the error in /2 derived in our calibration (see § 5 and Table 1). The calibration constant 
seems to be precise enough and does not limit the wavelength calibration. Finally, the solid 
line corresponds to a center uncertainty given by 6xc cos = 3 pixels. It produces an error in 
the wavelength scale that is almost 2 A in the outer parts of the FOV, exceeding the thresh- 
old set by the central wavelength. Three pixels represent a significant overestimate of the 
uncertainties we expect for the centers; if we consider a more realistic 6xc cos0 < 1.5 pixels, 
then 5\{R) < lA, and the error in the center position does not limit the calibration. We 
use this overestimated value to illustrate that the precision in the center turns out to be the 
critical point of the wavelength calibration. 

The question arises as to whether 5\{R) < lA suffices to compute the line fluxes that 
triggered our work (§1). We address the problem assuming the TF bandpass to behave as 
expected (e.g.. Born & Wolf 1980; Beckers 1998), 



r 



1 + 



2iXs-XiR)y '' 



with o"A the FWHM of the TF transmission bandpass, and X{R) the wavelength calibrated 
above. Consider a spectral line of wavelength As and integrated flux F. Then, the flux 
measured in a given point of the focal plane R, Fm{R), turns out to be 

Fm{R)=TF, (19) 

where we have assumed the line to be much narrower than the TF transmission (as it is 

indeed the case for the emission lines in typical HII regions). Equation (19) provides the 

recipe to retrieve the true flux from the observed one knowing the transmission. If X{R) is 

uncertain then r is uncertain as well, producing the error on F that we try to evaluate. The 

error in the retrieved flux SF due to the error in the central wavelength 6X{R) was derived 

by applying the law of propagation of errors (e.g., Martin 1971) to Eq. (19), and it is given 

by 

SF dlnr ^, , „, ,„„, 
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Fig. 7. — Wavelength errors in the calibration {5X) as a function of the distance to the 
center of the FOV. As indicated in the inset, the curves represent a mis-tuning of the central 
wavelength (5A(0) ^ 0, dotted line), errors in calibration constant {5C* ^ 0, dashed line), 
and uncertainties locating the center of the FOV {5xc 7^ 0, solid line). The errors we 
consider are realistic, except perhaps for an overestimate of 5xc (see main text). The central 
wavelength was assumed to be tuned to Ha. Vertical dotted lines show the radius of the MF 
(~2') and the non- vignetted FOV (~4'). 



After some trivial manipulations, the previous equation becomes 



Since < r < 1, then [(1 — r) r 



5F 
F 


~ -4[(1 - 


r)r] 


1/2 


6\{R) 


rl^/^ 


< 0.5, and 


so 
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6F 
F 


<2 


5\{R) 


(y\ 



(22) 



For 5\{R) < lA and a TF FWHM of a^ ~ 20 A, Eq. (22) yields \5F/F\ < 0.1. This 
error is small enough for most applications, including those mentioned in § 1 - e.g., it 
provides SFRs within 10% using the customary recipes (e.g., Kennicutt 1998), and it is 
much smaller than the intrinsic scatter of the sulphur doublet ([SII]AA6717, 6731) used to 
determine electron densities (e.g., Bresolin & Kennicutt 2002). Three final comments are 



-19- 



in order. First, a scaling factor r does not affect our error estimate, wliicli is tlie reason 
why the peak transmission in Eq. (18) was set to one. Second, Eq. (19) imphcitly assumes 
the observed spectrum to have no continuum. If there is continuum, then the arguments 
remain vahd replacing Fm{R) with the difference between the fluxes measured in line and 
in continuum. Finally, the upper limit in Eq. (22) significantly overestimates the true error 
when the wavelength setting is close to the peak transmission (see Eq. (21) with r ~ 1). 



6. Wavelength calibration protocol 

Based in our previous experience, this section puts forward a set of guidelines to carry 
out an accurate wavelength calibration of the upcoming scientific data. It involves routine 
observations of calibration lamps together with the astronomical images. Since no extensive 
tests on the stability of the OSIRIS TF center and calibration constant have been performed 
yet, they are needed to ensure the accuracy presented in this work and to identify other 
possible errors. The protocol will also allow to study and quantify the wavelength dependence 
of the calibration parameters. 

In general, the calibration process should consist of the following four steps: 

1. Obtain a number of calibration lamp images during the observation run covering the full 
FOV. Preferably, select a monochromatic emission line from a lamp with a wavelength 
near to the scientific observation. Then, map its positions along the OSIRIS FOV by 
changing the central wavelength of the TF. 

2. Calculate the center of the system using the rings corresponding to the monochromatic 
lines. Check that it is the same for all rings within a reasonable error (a few pixels; 
see Fig. 4). 

3. Calculate the value of the effective focal length of the camera lens (/2) by fitting linearly 
Eq. (12). This equation should be fitted in pixels or in a physical length (mm) to avoid 
crosstalk with errors in the plate scale. 

4. Calibrate the scientific images by assigning the wavelength given by Eq. (4) to every 
pixel. Use the empirically derived /2 and optical center. 

We recommend to apply the calibration described here to each CCD independently, 
and before applying astrometric corrections to avoid other error sources. Eventually, if the 
instrument turns out to be stable, the calibration will consist only of step 4 with /2 and the 
center derived in this paper. 
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For observations already in hand using the full OSIRIS FOV, we suggest to recalibrate 
the images using the prescriptions and constants given in this paper (Eqs. (4), (11), and 
(13)). Whenever a calibration ring is available as ancillary data, we recommend computing 
the center and compare its value with our estimate (Eq. (11)). A remainder is important: 
as we explained in §5.1, the trials for determining the best center have been carried out on 
the raw data. If the data to be re-calibrated have already been processed (e.g., overscan 
removed), then the values of the best center have to be shifted accordingly. 

Existing observations using only OSIRIS MP do not need re- calibration. The approxi- 
mate original prescriptions in OSIRIS manual suffice. 



7. Conclusions 

OSIRIS TF and GTC represent a unique combination to study the physical properties 
of extended astronomical sources. They are capable of providing high S/N, narrow band 
images, over a relatively large field with an excellent spatial sampling. The OSIRIS TF is 
based on a conventional FPI mounted in collimated beam. FPI based instruments provide 
an area of the FOV, the MP, where the wavelength can be considered constant within the 
band-pass of the TF. However, outside the MP there is a significant shift in wavelength which 
depends on the distance to the optical axis. Therefore, in order to benefit from OSIRIS large 
FOV, an accurate wavelength calibration is needed. 

Motivated by inconsistencies in the measurements of line ratios during an observational 
campaign, we decided to re-calibrate empirically the equation which describes the shift in 
wavelength with position on the FOV. Using spectral calibration lamps, we found that the 
proposed nominal calibration suffices only in the central 2'of the FOV, but fails elsewhere. 
The approximate equation originally given to describe the wavelength shift with radius is 
not correct for large radii, and the exact equation must be used instead (Eq. (4)). 

We identify three different sources of error in the wavelength calibration, namely, errors 
in the central wavelength (5A(0), errors in the calibration constant 6C*, and errors in the 
center position 6xc and 6yc. The error in the central wavelength is intrinsic to the procedure 
used to tune the instrument at the telescope. We conffim that the errors introduced in 
this process are < lA, as stated in the instrument manual. The calibration constant in the 
procedure proposed here only depends on the effective focal length of the camera lens. It 
has been empirically determined with a small error (Eq. (13)), good enough to discard the 
calibration constant as a source of uncertainty. We have also found a new optical center 
defined independently for the two OSIRIS CCDs (Eq. (11)) which improves significantly the 
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calibration. Actually, the position of the center seems to be the factor limiting the final 
accuracy. Considering all uncertainties, the updated calibration provides an accuracy of 
AA < lA in the entire OSIRIS FOV. 

The use of the new calibration makes OSIRIS ideal for studies involving spectropho- 
tometric measurement of large galaxies covering the full FOV, such as those that triggered 
the present work (see § 1). However, it is of broader interest. Whatever study using the 
entire FOV will benefit from the new calibration (e.g., extended planetary nebulae, clusters 
of stars or galaxies distributed in a large area, etc). 

Based in our previous experience, a wavelength calibration protocol is proposed in § 6. 
It will allow to check the stability of both OSIRIS center and the calibration constant, i.e., 
to see whether they vary over long periods, when OSIRIS optics is re-aligned, or if different 
wavelength ranges are used. In particular, when the blue etalon will be at work, a similar 
study on the wavelength calibration must be done. 

Thanks are due to the GTC staff, in particular to A. Cabrera-Lavers, for their support 
and for providing the data used in this study. Discussions with J. L. Rasilla, R. Lopez, 
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ICM/GTC, were extremely clarifying. We are also indebted to A. Manescau and J. A. Bonct 
for support and advice. This work has been partially funded by the Spanish MICINN (grants 
AYA2007-67965-C03-01 and AYA2010-21887-C04-04) and the Consohder-Ingenio Program 
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